The N to A electromagnetic transition form factors from Lattice QCD 
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The magnetic dipole, the electric quadrupole and the Coulomb quadrupole amplitudes for the 
transition 7A'' — > A are calculated in quenched lattice QCD at /3 = 6.0 with Wilson fermions. 
Using a new method combining an optimal combination of interpolating fields for the A and an 
overconstrained analysis, we obtain statistically accurate results for the dipole form factor and for 
the ratios of the electric and Coulomb quadrupole amplitudes to the magnetic dipole amplitude, 
Rem and Rsm, up to momentum transfer squared 1.5 GeV^. We show for the first time using 
lattice QCD that both Rem and Rsm are non-zero and negative, in qualitative agreement with 
experiment and indicating the presence of deformation in the N- A system. 
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Deformation is an important and well studied phe- 
nomenon in atomic and nuclear physics, and it is de- 
sirable to understand whether it also arises in low-lying 
hadrons and if so, why. For classical and quantum sys- 
tems with spins larger than 1 /2, the one-body quadrupole 
operator provides a convenient characterization of defor- 
mation. Experimentally, however, in the excited spin 3/2 
A, which can have a non-zero quadrupole moment, it is 
not practical to measure it, and the spin 1/2 nucleon, 
which is easily accessible to measurement, cannot have 
a spectroscopic quadrupole moment. Hence, the experi- 
ment of choice to reveal the presence of deformation in 
the low-lying baryons is measuring the N -A transition 
amplitude, and significant effort has been devoted to pho- 
toproduction experiments on the nucleon at Bates and 
Jefferson Lab in order to measure to high accuracy the 
ratios of the electric (E2) and Coulomb (C2) quadrupole 
amplitudes to the magnetic dipole (Ml) amplitude. If 
both the nucleon and the A are spherical, then E2 and C2 
are expected to be zero. Although Ml is indeed the domi- 
nant amplitude, there is mounting experimental evidence 
over a range of momentum transfers that E2 and C2 are 
non-zero [3j- Similarly in lattice QCD, for hadrons with 
spins larger than 1/2, the deformation is determined by 
measuring their quadrupole moment knowing the hadron 
wave function, which can be obtained via density correla- 
tors Using these techniques, it was shown that the 
rho has a non-spherical spatial distribution with a non- 
zero quadrupole moment and that the A acquires a small 
deformation as the quark mass decreases However, 
direct contact with experiment is established by calcu- 
lating the N to A transition form factors. 

In this work we calculate these form factors as a func- 
tion of the momentum transfer in lattice QCD in the 
quenched approximation on a lattice of size 32"^ x 64 at 
(3 = 6.0. We obtain, for the first time, accurate results for 
the E2 and C2 moments for momentum transfer squared, 
g^, up to about 1.5 GeV^. Our results are sufficiently ac- 



curate to exclude a zero value. The two novel aspects, as 
applied to the N to A matrix elements, that are crucial 
for obtaining this accuracy are: 1) An optimal combi- 
nation of three-point functions, which allows momentum 
transfers in a spatially symmetric manner obtained by an 
appropriate choice of the interpolating field for the A. 2) 
An overconstrained analysis using all lattice momentum 
vectors contributing to a given value in the extraction 
of the three transition form factors 

In lattice QCD, transitions involving one-photon ex- 
change, such as the one involved in the N to A transi- 
tion, require the evaluation of three-point functions. The 
standard procedure to evaluate a three-point function is 
to compute the sequential propagator. This can be done 
in two ways: In the first approach used in previous lat- 
tice calculations 0)0) the photon couples to a quark at a 
fixed time ti carrying a fixed momentum q. This means 
that the form factors can only be evaluated at one value 
of the momentum transfer. Since the current must have 
a fixed direction and a fixed momentum this approach 
is referred to as the fixed current approach. Within this 
approach one can use any initial and final state without 
requiring further inversions, which are the time consum- 
ing part of the evaluation of three-point functions. In 
the second approach, which is used in this work, we re- 
quire that the initial state created at time zero has the 
nucleon quantum numbers and the final state, annihi- 
lated at a fixed time t2, has the A quantum numbers. 
The current can couple to any time slice ti carrying any 
possible value of the lattice momentum |^ |^ . Because 
the quantum numbers of the final state are fixed we re- 
fer to the second method as the fixed sink method. With 
the improvements implemented in this work, this method 
becomes clearly superior to the fixed current approach al- 
lowing accurate evaluation of the form factors as a func- 
tion of q^. 

The matrix element for the — > A transition with 
on-shell nucleon and A states and real or virtual photons 
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has the form (T3| 

(A(p',s') |j^|7V(p,s)) = 

1/2 



/2 
iM - 



urip',s')0^^u{p,s) (1) 



where p(s) and p'{s') denote initial and final momenta 
(spins) and Ur{p',s') is a spin- vector in the Rarita- 
Schwinger formalism. The operator can be decom- 
posed in terms of the Sachs form factors as 

^ g,niq')Kl^,+gE2iq')KZ + gc2iq'Wc2 - (2) 

where the magnetic dipole, Gmi, the electric quadrupole, 
gE2, and the Coulomb quadrupole, t/c2, form factors de- 
pend on the momentum transfer = (p' — pY' ■ The 
kinematical functions K"^^ in Euclidean space a re g iven 
in ref. 8]. Using the relations given in refs. [lOl the 
ratios Rem and Rsm in the rest frame of the A are ob- 
tained from the Sachs form factors via 



Rem — — 



QE2{q^) 

GMi{q^) 



Rsm — — 



|q| Qc2{q^) 

2mA Gmi{,cP) 



(3) 



The ratio Rem is also known as EMR and Rsm as CMR. 

To extract the N to A matrix element from lattice 
measurements we calculate, besides the three point func- 
tion G^^"^ {t2,ti;p ', p; r), the nucleon and A two-point 
functions, and G^^ and look for a plateau in the 

large Euclidean time behavior of the ratio 

{G^^"^(t2,ti;p',p;r)> 



1/2 



(4) 



(G^^(t2 - ti, p; r4)) (G^^(ti, p r4)) (G,^^(t2, p '; r4)) 
(G^^(f2-ti,p';r4)) (G'^^(ii,p;r4)) (G~^^(t2,p;r4)) 

^ n<j(p ,p ;r;^) . 

We use the lattice conserved electromagnetic current, 
j^{x), symmetrized on site x by taking j^{x) — > 
+ — jl)] /2 and projection matrices for the 
Dirac indices 
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Throughout this work we use kinematics where the A 
is produced at rest and therefore q = p' — p = —p. 
We fix i2 = 12 in lattice units and search for a plateau 
of i?a-(i2, ^i; P P as a function of ti. — —q^ 

denotes the Euclidean momentum transfer squared. 

We can extract the three Sachs form factors from the 
ratio of Eq. by choosing appropriate combinations 
of a indices and F matrices. However there are several 
choices of a indices and F matrices that can be used each 
requiring an inversion. Therefore we must determine the 
most suitable combination of three-point functions from 
which to extract the Sachs form factors. 

For example the dipole form factor can be extracted 
from 



n,(q ;r4;/i) = lAe^^^^p^GMiiQ^) 



where A is a kinematical coefficient. This means that 
there are six statistically independent matrix elements to 
extract Gmi each requiring the evaluation of a sequential 
propagator. However, due to the epsilon factor, a choice 
of one of the six combinations means that only momen- 
tum transfers in one direction contribute. Instead if we 
take the symmetric combination. 



Siiciifi) = ^n^(q ;F4;^) 



(7) 



lattice momentum vectors in all directions contribute. 
This combination, which we refer to as sink type Si, is 
built into the A interpolating field and requires only one 
inversion. 

Another choice of three-point functions is to use the 
projection matrices F^ instead of F4. The relations are 
more involved in this case and will be discussed in detail 
in a forthcoming publication [l^ . However, as in the 
example given in Eq. ©, for any current direction, one 
can choose any of the six statistically different three-point 
functions from which to extract the Sachs form factors. 
Instead of choosing one of six we can consider a linear 
combination that involves, in a symmetric manner, all 
spatial directions allowing, for a given Q^, the maximum 
number of momentum vectors to contribute. We take 



S2ici; fJ.) = n^(q iTfei^) 



(8) 



which we refer to as sink 82- When the current is 
in the spatial direction both Ge2 and Gc2 can be ex- 
tracted from 5*2 with one inversion. In addition, when 
the current is in the time direction, S2 provides a statis- 
tically independent way for evaluating Gc2, with no ex- 
tra cost. Another combination to extract Ge2 and Gc2 is 
^3 (q; = Ha (q; F3 ; m) - (Hi (q; Fi ; /i) + H2 (q; F2 ; m)) /2, 
which produces results of comparable quality to those ob- 
tained with 5*2 13]. In the case of E2, sink type S3 has 
the disadvantage of vanishing at the lowest value of 
whereas 52 contributes at all values of Q^. For C2, on 
the other hand, 6*2 gives zero at the lowest Q^, whereas 
source type S3 gives a non- vanishing result. 

The second important ingredient in the extraction of 
the form factors is to take into account in our analysis 
all the lattice momentum vectors that contribute to a 
given Q^. This is done by solving the overcomplete set 
of equations 



P{cl;^l) = D{q■,^i)■F{Q^) 



(9) 



(6) 



where P(q; /i) are the lattice measurements of the ratio 
given in Eq. Q having statistical errors Wk and using the 

(Gmi \ 
Ge2 and, with N being 
GC2 J 

the number of current directions and momentum vectors 
contributing to a given Q^, D is an A'^ x 3 matrix which 
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depends on kinematical factors. We extract the form 
factors by minimizing 




(10) ^ 



using the singular value decomposition of D. 
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FIG. 1: Gmi as function of Q^. Results a.t k — 0.1554 are 
shown by the crosses, at k = 0.1558 by the open triangles 
and at K = 0.1562 by the asterisks. The filled circles are 
results in the chiral limit and the filled triangles are results 
extracted from fits to the measured differential cross sections 
usmg MAID [a. The filled square is the Particle Data Group 
result flHI . The solid lines are fits using the Ansatz of Eq. II12|I . 
The dashed line is a fit to the lattice data using the Ansatz 
aexp(— feQ^). 



All the results for the form factors are obtained us- 
ing 200 configurations and three values of the hopping 
parameter k. The values of k chosen are 0.1554,0.1558 
and 0.1562 and give ratio of pion to rho mass rriT^/mp — 
0.64, 0.59 and 0.50 respectively. We use the nucleon mass 
at the chiral limit to set the lattice spacing a obtaining 
= 2.04(2) GeV (a=0.098 fm). Using the optimized 
sink Si we show in Fig. ^ at the three quark masses our 
results for the Ash form factor Q^n defined by [l^ 
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(mjv+mAp' 



(11) 



To obtain the results at the chiral limit, shown on the 
same figure, we perform a linear extrapolation in m^. 
Although we expect chiral logs that appear at next-to- 
leading order in chiral perturbation theory to be sup- 
pressed for the momentum transfers studied in this work, 
our linear extrapolation introduces a systematic uncer- 
tainty. This uncertainty can not be assessed, since the 
known chiral perturbation theory results jlllll7l | are valid 
at very low masses and momentum transfers. On the 
same figure we also show the experimental values as ex- 
tracted from the measured cross sections using the phe- 
nomenological model MAID . We perform fits to both 
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FIG. 2: Rem as a function of . Upper graph shows re- 
sults at «: = 0.1554 (crosses), at k = 0.1558 (open triangles) 
and at K = 0.1562 ( asterisks). The filled circles show chiral 
limit extrapolations. Lower graph shows The filled triangles 
show recent experimental results |l| (filled triangles) and Q 
)opensqures). The filled square is the Particle Data group 
result [il. 

the lattice data at the chiral limit and to the MAID data 
using the phenomenological parametrization 



^a(Q') =ea(0)i?a(Q')G|(Q') 



(12) 



where Ra{Q^) for a — Ml, E2 and C2 measures the de- 
viations from the proton electric form factor G^^{Q'^) = 
1/(1 + Q^/0.71)^. Usually experimental data are fit- 
ted by taking RjpiQ^) = Re2{Q^) = Rc2{Q^) = 
l-|-aexp(— 7(5^) [l8i|. As can be seen this fitting Ansatz 
provides a very good description to the MAID data. Al- 
though the lattice data at the chiral limit lie higher than 
the MAID data they can be fitted to the same form yield- 
ing at = a value consistent with that given by the 
Particle Data Group j^. The lattice data are also well 
described by the simple exponential Ansatz a exp(— 6Q^), 
which however at = gives a value lower than exper- 
iment. 

Using the optimized sink 6*2 we extract the quadrupole 
form factors, Ge2 and Gc2, at three values of the quark 
mass, as shown in Figs. |21 and |3| We note that Qc2 at 
the lowest value is extracted using source type S3 
since S2 vanishes for this particular lattice vector. On 
the same figures we also show the values obtained in the 
chiral limit by performing a linear extrapolation in mj. 
As expected, both EMR and CMR become more neg- 
ative as we approach the chiral limit. Our results for 
EMR and CMR at the chiral limit are compared to recent 
measurements [3, 01 in Figs. |21 and |21 respectively. The 
quenched results for EMR are accurate enough to exclude 
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FIG. 3: i?sM as a function of Q^. The notation is the same 
as that of Fig. H 

a zero value at low Q^. With our current statistics they 
are in agreement with the experimental measurements. 
Whether the apparent discrepancy for CMR at low 
is a significant deficiency of quenched QCD or a problem 
with a single data point remains to be resolved by new 
measurements that are currently being analyzed. As 
increases, however, the quenched results agree qualita- 
tively with measurements. 

In summary, two novel methods are applied in the eval- 



uation of the N to A transition form factors: The first 
improvement comes from employing an optimized sink 
for the A allowing a maximum number of lattice ma- 
trix elements to contribute and the second from utilizing 
this enlarged set of data in an overconstrained analy- 
sis. Given that there are ambiguities in the extraction 
of the quadrupole amplitudes from experimentally mea- 
sured response functions arising from using models, accu- 
rate lattice data are extremely valuable: excluding a zero 
quadrupole stength in lattice QCD corroborates experi- 
mental observations for a non-zero Rem and Rsm- This 
is particularly important for CMR where at low there 
are very few accurate experimental measurements. If 
confirmed, the agreement of EMR with experiment at low 
while CMR disagrees with experiment raises interest- 
ing questions regarding the pion cloud contributions to 
these ratios due to the absence of the sea quarks. Hav- 
ing, for the first time, demonstrated that we can reliably 
extract CMR in quenched lattice QCD opens the way 
for a precise investigation of sea quark contributions to 
both EMR and CMR, which can lead to an understand- 
ing of the physical mechanism responsible for non-zero 
quadrupole strength in the N to A transition. 
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